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1. Introduction 

In the last years there has been a growing interest in nonequilibrium phase transitions 
[T|,|2]. The absence of a general theory nonequilibrium systems originates a great number 
of open problems, even in one-dimensional systems that, in the equilibrium regime, 
generally are exactly solvable. Usually, numerical simulations are a useful technique 
in the study of phase transitions and critical phenomena, and its power has been 
growing with the increasing capacity of the computers and the development of new 
simulational techniques. Non-equilibrium models are particularly suited for simulations. 
However, others approaches may be complementary in the study of these phenomena. 
Thus, it is interesting study the systems by other techniques. Among these alternative 
methods we can cite the series expansions [3], which in some cases lead to very precise 
estimates for the critical properties that characterize these transitions [3], H] and the 
exact diagonalization of the time evolution operator in finite-size systems [5]. 

A special case of the non-equilibrium models are the ones that exhibit absorbing 
states [1], that is, states that may be reached in their dynamics, but transitions 
leaving them are forbidden. Since such models do not obey detailed balance, they 
are intrinsically out of equilibrium. The most studied system for this class of problems 
is the so called contact process (CP) model [6], a toy model for the spreading of an 
epidemic. This model displays a transition between an absorbing and an active state 
with critical exponents belonging to the directed percolation (DP) universality class [7]. 
In addition, the CP model is related to the Schlogl's lattice model for autocatalytic 
reactions [8J and the Reggeon Field Theory [9J. 

Many variants of this model have been studied flU{ [TT| [T2| [13] . most of them 
belonging to DP class also. In fact, the robustness of this universality class is an evident 
characteristic of these models. Such robustness is explained by the conjecture that all 
models with phase transitions between active and absorbing states with a scalar order 
parameter, short range interactions and no conservation laws belong to this class [Ti] . 

One of these variants is the CP with diffusion |15], which exhibits a critical line 
instead of a critical point. This line begins at the critical point of the model without 
diffusion and ends in the infinite diffusion rate limit, where the critical properties of the 
system approaches those predicted by the mean-field approximation. For finite values of 
the diffusion rate, the critical behavior of this model is dominated by the DP universality 
class. This is not surprising, since the original dynamic includes an intrinsic diffusion 
process. On the other hand, the mean-field behavior in the limit of infinite diffusion may 
be understood considering that, since diffusion processes are dominating the evolution 
of the system in this limit, creation processes are effectively determined by the mean 
densities, as is supposed in the mean-field approximation. This change of behavior at 
the infinite diffusion rate limit, between the critical behavior of the DP class and the 
one predicted by the mean-field approximation, characterizes a crossover of the critical 
properties in the neighborhood of a multicritical point. As in the equilibrium case we 
may then write any density variable, in the neighborhood of a multicritical point, as the 
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following scaling form [TB]: 



g{a,D) = {a.-afF 



Dc-D 



(1) 



where a is a transition rate for annihilation of particles, D is the diffusion rate, 
ranging between and 1, ^ is a critical exponent associated to the density variable 
(yf, corresponding to the value predicted by the mean-field approximation and is the 
crossover exponent. The scaling function F{z) is singular at a point z = zq oi its 
argument, such that the critical line, in the neighborhood of the multicritical point, 
corresponds to 



One of the first studies of this problem was performed by Dickman and Jensen 
[T5] . who considered the model using supercritical series in a with the diffusion rate 
D taken as a fixed parameter. Therefore, in their calculation series expansions are 
derived for fixed values of the diffusion rate and the analysis of this series leads, 
among other information, to the phase diagram of the model with diffusion. However, 
they found that the fluctuations of the estimates provided by d-log Pade approximants 
increase as the diffusion rate grows, so that the critical curve was obtained only up to 
D 0.8. Since the crossover exponent characterizes the critical curve close to the 
infinite diffusion rate limit D ^ 1 no precise estimate of the crossover exponent was 
possible. The disappointing performance of the Pade approximants as the multicritical 
point is approached is not surprising, since it is known tha one-variable series analysis 
techniques are not effective close to such points [16j. More recently [IT], the model was 
simulated in a conservative ensemble. These simulations display smaller fluctuations in 
the estimates, making it possible to obtain the critical line up to values close to the 
multicritical point, furnishing the value (f) = 4.03(3) for the crossover exponent. This 
result is consistent with the lower bound > 1 predicted by Katori in |17j . 

In the present work, we obtain estimates for the critical line, exponents (3 and z 
as well as the crossover exponent for the one-dimensional CP with diffusion. For this 
task we use two approaches: a two-variable supercritical series and exact solutions 
for flnite-size one-dimensional lattices. The supercritical series is analyzed using a 
partial differential approximants (PDAs) flGl fT8| 119]. This technique seems to be more 
appropriate for the analysis of a two-variable series with a multicritical behavior, as is 
shown by the results obtained for other models [121 EQ]- This paper is organized as 
follows. In section II we present the model and the mean-fleld results, in section III we 
show the derivation of the supercritical series and in section IV the analysis of this series 
is presented. In section V we discuss the flnite-size exact solution and their results for 
diffusive CP model. Finally, in section VI, the conclusions and flnal discussions of this 
work may be found. 




(2) 
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2. Definition of the model and mean-field results 

In a one-dimensional lattice each site can be empty or occupied by a particle, so that 
we will associate an occupation variable r/j = 0, 1 to the site i. The evolution of the 
system is governed by markovian local rules such that the particles are annihilated with 
rate a and created in empty sites with a transition rate n/2, where n is the number of 
occupied nearest neighbors and z is the total number of nearest neighbors. In addition 
to these rules that define the CP, we include a diffusive process, allowing the hopping of 
particles to empty nearest neighbor sites at the rate D = D/{1 — D). The configuration 
such that all sites are empty is an absorbing state. The passage from an active steady 
state, with a nonzero density of particles, to the absorbing state defines a transition line 
in the {a, D) plane as shown in figure [H 

A mean-field approach for this model can be obtained at several levels of 
approximation [Ij. In the one-site level the role of the diffusion is irrelevant since it 
contributed equally to creation and annihilation of particles at a given site i. Already 
in the two-site level it is possible to determine the critical line by using as variables the 
parameters a and D. This line is described by the expression 



and is very easy to show that in the neighborhood of the multicritical point, = 
1,D = 1), the behavior of this curve is given by the scaling relation 



where = 1. On the three-site level an unitary crossover exponent also appears, as may 
be seen in figure [H 

This result is in accordance with the lower bound, > 1, determined by Katori |17j . 
Comparing the mean- field approximation with the simulational result [11], we observe 
that the first approach always overestimates the supercritical region of the models. Even 
in the higher diffusion region, the critical line obtained using mean-field does not coincide 
with the simulational result. Actually, this coincidence occurs only in the multicritical 
point. 

In the next section we will derive the supercritical series in the variables A and D 
to determine the value of this exponent and compare it with that obtained in [TT], [TB] . 

3. Derivation of the supercritical series for the model 

We use the operator formalism proposed by Dickman and Jensen [1] in order to derive a 
supercritical series. To this end, we define the microscopic configuration of the lattice, 
|?7), as the direct product of kets |?7) = 0j \r]i), with the following orthonormality 
property, {r]\r]') = YliSff^^rn/- The particle creation and annihilation operators at site 
i are defined as 




(3) 




(4) 



(5) 
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Figure 1. Top: phase diagram obtained using the two- and three-site mean- field 
approach and results of simulations [11] . Bottom: the log- log plot of the same quantity 
is plotted and compared to the result (f) — 1 and (/> = 4 for these approaches. 



In this formalism, the state of the system at time t is 

\m) = j:piv,t)\v), (6) 

where p{ri, t) is the probabihty of a configuration rj at time t. If we define the projection 
onto all possible states as ( | = J2{ri}{v\ then the normalization of the state of the system 
may be expressed as ( {ip) = 1. In this notation, the master equation for the evolution 
of the state is: 

*^ = S|V,W>. (7) 
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The evolution operator S may be expressed in terms of the creation and annihilation 
operators as S' = fiR + V where 

R = DY.il- 4_iA)v4i_i4 + (1 - 4+iA)A+i4 + 

i 

+ Y^{A,-AlA,), (8) 

i 

V = Y^iAl - aaI){aUa_^ + 4+i^.+i), (9) 

i 

where /i = 2a and D = D/a. 

We notice that the operator R diffuses (01 — > 10) or annihilates particles (1 0), 
while the operator V acts in the opposite way, generating particles (0 — > 1). It is 
convenient to join the diffusion with the annihilation process to avoid ambiguities in 
the truncation of the series at a certain order. For small values of the parameter /i 
the creation of particles is favored, and the decomposition above is convenient for a 
supercritical perturbation expansion. Using the equations ([8]) and ([9]) the action of each 
operator on a generical configuration (C) is given by 

R{c) = D E(c:) + E(cr) + E(ci') + 

_ i j ^ . 

+ E(Cn-[(n + 2r2) + D + r](C), (10) 

i 

where the first sum is over ri sites with particles and one empty neighbor, the two next 
sums are over r2 sites with particles and two empty neighbors and the last sum is over 
all sites occupied by a particle. Configuration (C^') is obtained moving the particle at 
the site i to the single empty neighbor site, (C-^*^''^) is a configuration where the particle 
at the site i moved to the empty neighbor at the right (r) or at the left (/) is replaced 
by a hole and one of the empty neighbors (at the right or left) is occupied. On the other 
hand, the action of operator V is 

V{C) = YSC'n + 2E(C;') - (gi + 2q,){Cl (11) 

« 3 

where the first sum is over the qi empty sites with one occupied neighbor, the second 
sum is over the q2 empty sites with two occupied neighbors. Configuration (Cf) is 
obtained occupying the site i in configuration (C). 

To obtain a supercritical expansion for the ultimate survival probability of particles, 
we start by remembering that in order to access the long-time behavior of a quantity, 
it is useful to consider its Laplace transform, 

ms)) = / e-^'m))- (12) 
Jo 

Inserting the formal solution li'it)) = e'^*|V'(0)) of the master equation (j7]) we find 

ms)) = {s-S)-'\m)- (13) 
The stationary state \ip{(X))) = lim^^oo \'ip{t)) niay then be found by means of the relation 
|V^(oo)) = lims|Vi(s)). (14) 
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A perturbative expansion may be obtained by assuming that {ipis)) can be expressed in 
powers of /i and using f|T3l) . 

ms)) = i^po) + + ^'\^2) + ... = {s-v- iiR)-'\m)- (15) 

Since 

{s-V-jjR)-^ = (s-V)-^ [l + - Vy^R + - Vy'^R'^ + ■■■], (16) 

we arrive at 

lV;o) = (^-V^)-^|^(0)) (17) 

and 

\^n) = {s-V)-'R\tPn-l), (18) 

for n > 1. The action of the operator (s — V)~^ on an arbitrary configuration (C) may 
be found by noticing that 

(. - Vr\C) = s~' {(C) + (. - V)-'V{C)} , (19) 

and using the expression [11] for the action of the operator V, we get 



s-V)-\C) = sAiC) + is-V) 



E(cr)+2E(cr) 



(20) 



where the first sum is over the qi empty sites and one occupied neighbor, the second sum 
is over the q2 empty sites and two occupied neighbors, and we define Sg = l/(s+gi + 2g2)- 

It is convenient to adopt as the initial configuration a translational invariant one 
with a single particle (periodic boundary conditions are used). Now, looking at the 
recursive expression ( !20l) . we may notice that the operator [s — V)~^ acting on any 
configuration generates an infinite set of configurations, and thus we are unable to 
calculate lip) in a closed form. However, it is possible calculate the extinction probability 
p{s), which corresponds to the coefficient of the vacuum state |0). As happens also for 
models [H [3] related to the CP, configurations with more than j particles only contribute 
at orders higher than j, and since we are interested in the ultimate survival probability 
for particles Poo = 1 — lim^^o -5^(5), Sg may be replaced by 1/q in equation fl20l) . An 
illustration of this procedure may be found in a previous calculation [12]. 

The algebraic operations described above is performed by a simple algorithm which 
allow us to calculate 24 terms with a processing time to about 2 hours. Actually, the 
limiting factor in this operation is the memory required. In this way we define the 
coefficients bi j for the ultimate survival probability as: 

11 24 i-2 

Poo = 1 - - -/.^ - ^ 5: k,i/D\ (21) 

^ ^ i=3 j=0 

and they are listed in Table [TJ 
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Table 1. Coefficients for tlie series expansion for ultimate survival probability 
corresponding to the CP model with diffusion. The indexes refer them to the equation 



4 



6 



7 




1 


1 

2 

1 
2 
3 

1 
2 
3 
4 

1 
2 
3 
4 
5 

1 
2 
3 
4 
5 
6 

1 
2 
3 
4 
5 
6 
7 

10 
1 
2 
3 
4 
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0.25000000000000000000 x 10° 
-0.25000000000000000000 x 10° 
0.28125000000000000000 x 10° 
-0.37500000000000000000 x 10° 
0.37500000000000000000 x 10° 
0.34375000000000000000 x 10° 
-0.50781250000000000000x10° 
0.57812500000000000000x10° 
-0.62500000000000000000 x 10° 
0.44726562500000000000 x 10° 
-0.76220703125000000000x10° 
0.85058593750000000000 x 10° 
-0.83984375000000000000 x 10° 
0.10937500000000000000x10° 
0.60223388671874955591 x 10° 
-0.11734619140625004441 xlO^ 
0.15190429687499997780x10^ 
-0.14140624999999984457x10^ 
0.10878906249999982236x10^ 
-0. 19687500000000000000 x 10^ 
0.83485031127929687500x10° 
-0.18110389709472716202x10^ 
0.25234603881835981909 x 10^ 
-0.29291381835937464473 x 10^ 
0.24864501953125062172x10^ 
-0.10754394531250017764x10^ 
0.36093750000000000000 x 10^ 
0.11814667913648828623x10^ 
-0.28569926950666579835 x 10^ 
0.42781094621729156557x10^ 
-0.48761836864330092567x10^ 
0.54410674483687788694 xlO^ 
-0.48496839735243062464 xlO^ 
0.11848958333333414750x10° 
-0.67031250000000000000 xlO^ 
0. 16988672076919952847X 10^ 
-0.45030223008843064392 x 10^ 
0.73700355965291270977x10^ 
-0.92486491500105252328 x 10^ 
0.87502182305104447835 x 10^ 



11 



12 



13 



6 
7 
8 

1 
2 
3 
4 
5 
6 
7 



14 




1 
2 
3 
4 
5 
6 
7 
8 
9 

10 



1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 



1 

2 



0.10567643059624561630x102 
0.34998474121093847700X 10^ 
0.12568359375000000000X 102 
0.24775957118666096513x10^ 
-0.71703082845177412707X 10^ 
0.12429764028158224676x102 
-0.16676337106809967281 X 10^ 
0.19435605471721252968x102 
-0.15230894658300556443X 102 
0.13923689787279895924X 102 
-0.24910518081099844778X 102 
-0. 13853664539478481643 X 102 
-0.23740234375000000000 X 102 

0.36488812342264926869x10^ 
-0.11443929729648042226X 102 
0.21418735868689896762X 102 
-0.29831307350977681381 X 102 
0.34272964785342651339X 102 
-0.40398672142671166796X 102 
0.27855684964608855125X 102 
-0.13595902518316915319X 102 
0.60935236387946176251 x 102 
0.40270442479922508028 X 102 
0.45106445312500000000X 102 
0.54293656084851154020x10^ 
-0.18322144692814863021 x 102 
0.36259195896082665911x102 
-0.54866931326313050477X 102 
0.67130818799164941879X 102 
-0.64293858561553619779X 102 
0.81638245065997452343X 102 
-0.59206203158103356543X 102 
-0.11386872839957611347x102 
-0.15017302299755911577X 10^ 
-0.10358591484729181786x10=* 
-0.8611230468750000000X 102 
0.8132542219307161701600x10^ 
-0.29467694610727896531 x 102 
0.62075441392908530247X 102 
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i 3 hi 

14 4 0.12740413805065173847x10=^ 

5 -0.14872799806680012580x10^ 

6 0.11055259565281477308x103 

7 -0.15322889380857196784x103 

8 0. 15125687744198603468 xlO^ 

9 0. 12100775973170790678 xlO^ 

10 0.36717258144235222517x10^ 

11 0.24952145042880511028x10^ 

12 0. 16504858398437500000 xlO^ 

15 0.12275012836144505002 xlO^ 

1 -0.47363165128788978109x102 

2 0.10546586796137503939x103 

3 -0.1762739398241103288x103 

4 0.23298609118631188153x103 

5 -0.27071715385838172097x103 

6 0.33267266081610591755x103 

7 -0.18014432368848466126x103 

8 0.24323790718884771422x103 

9 -0.43351480900655008099x103 

10 -0.48110863117407200207x103 

11 -0.88513499744041246231x103 

12 -0.57707323452079549497x103 

13 -0.31740112304687500000x103 

16 0.18620961415130427241 xlO^ 

1 -0.76547748518027589171 xlO^ 

2 0.17936794520034777634x103 

3 -0.30967915791812674797x103 

4 0.45127497217248452444x103 

5 -0.53698491310493750461x103 

6 0.51912309945306844838x103 

7 -0.74776029427868388666x103 

8 0.31390176074091152714x103 

9 -0.23066668763466492464 x 103 

10 0.12806582574012411442x10^ 

11 0.15251186818533849419x10^ 

12 0.21006816788719602300x10^ 

13 0.12983815231244370807x10^ 

14 0.61213073730468750000x103 

17 0.28405733950686048672x10^ 

1 -0.12342415559750365617x103 

2 0.30541526863109334045x103 



3 -0.55737662171810006839x103 

4 0.80584299514354984240x103 

5 -0. 10663609887965685630 xlO^ 

6 0.12359653482806177180x10^ 

7 -0.86967640846259655518x103 

8 0.16397080915531578285x10^ 

9 -0.72176857313912660175x103 

10 -0.39467789553376610456x103 

11 -0.37241724510229601037x10^ 

12 -0.43433579993156563432x10^ 

13 -0.49111545112523144780x10^ 

14 -0.28643360597728060384x10^ 

15 -0.11834527587890625000x10^ 

18 0.435207828742268674200 xlO^ 

1 -0.20009555228747112210x103 

2 0.51666085550451919062x103 

3 -0.98221829678260564833x103 

4 0.15468813789798582548x10^ 

5 -0.19452358996876919264x10^ 

6 0.23142581042234874076x10^ 

7 -0.29269143532222028625x10^ 

8 0.11929249112512670763x10^ 

9 -0.33469260525748682085 xlO^ 

10 0.22662397929431922421x10^ 

11 0.33534053058169365613x10^ 

12 0.10497352179518215053x10^ 

13 0.11603782689493993530x10^ 

14 0.11326231062527707763x10^ 

15 0.62269615466431168898x10^ 

16 0.22929397201538085938x10^ 

19 0.66930218067969633466x10^ 

1 -0.32354897975245813768x103 

2 0.88042629554806353553 x 103 

3 -0.17413898514485581472x10^ 

4 0.27574673463423659996x10^ 

5 -0.39760107863198809355x10^ 

6 0.45478180647223589403x10^ 

7 -0.44829538123020120111x10^ 

8 0.71439950475520599866x10^ 

9 -0.11433220550468702186x10^ 

10 0.58699584417021997069x10^ 
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i 


j 




i 


; 


Kj 


19 


11 


-0.800073099504771 19855 x 10^ 




15 


-0.19550932821134978440x10^ 




12 


-0.13849586891129882133x10^ 




16 


-0.17865808405461237999x10^ 




13 


-0.28630836746692133602x10^ 




17 


-0.13001245697926016874x106 




14 


-0.29699369329023520550x10^ 




18 


-0.60334798749708410469 X 10^ 




15 


-0.25808082841629722679x10^ 




19 


-0.16853935146331787109x10^ 




16 


-0.13385541065918856475X 10^ 


22 





0.24876519640902955643x10^ 




17 


-0.44510006332397460938x10^ 




1 


-0.13849806980532957823x10^ 


20 





0. 1033739990888301 1790 X 10^ 




2 


0.42491525612070690840x10^ 




1 


-0.52548922262251915072 xlO^ 




3 


-0.95509490156905540061 x 10^ 




2 


0.14837268484319015442x10^ 




4 


0.17072122607398174296x10^ 




3 


-0.30794482965317188246x10^ 




5 


-0.25259991312968326383x10^ 




4 


0.51825134723219762236X 10^ 




6 


0.36959825267281092238x10^ 




5 


-0.70384028435684167562X 10^ 




7 


-0.39724946702482979163x10^ 




6 


0.95392163017937873519X 10^ 




8 


0.41082488201681495411x10^ 




7 


-0.10775090683332531626x10^ 




9 


-0.69996295461125846487X 10^ 




8 


0.72306238600702890835 x 10^ 




10 


-0.47433844045513687888X 10^ 




9 


-0.17556031650766155508x10^ 




11 


-0.94472082453477501986 X 10^ 




10 


0.45576565336850308086 X 10^ 




12 


0.18231631734823597071 x 10^ 




11 


-0.68853229355527937514X 10^ 




13 


0.74507403911089320900x10^ 




12 


0.27725591978402942914x10^ 




14 


0.28127919612702319864x10^ 




13 


0.46620453043310422800x10^ 




15 


0.40344708716105029453x10^ 




14 


0.75787894589888033806x10^ 




16 


0.49346542769156675786x10^ 




15 


0.73699439769879478263x10^ 




17 


0.42523410507562680868x10^ 




16 


0.58191226154708187096 x 10^ 




18 


0.28817317888963740552 x 10^ 




17 


0.28519783989193914749 xlO^ 




19 


0.12690253004534545471 x 10^ 




18 


0.86547234535217285156x10^ 




20 


0.32865173535346984863x10^ 


21 





0.15998830271612598608x10^ 


23 





0.38696067609131841891 x 10=^ 




1 


-0.85206464313625656359x10^ 




1 


-0.22521288931067451813x10^ 




2 


0.2527041 1365871 177622 x 10^ 




2 


0.72239403631839231821x10^ 




3 


-0.53883333002320678133X 10^ 




3 


-0.16540644823073282168x10^ 




4 


0.93069165541648508224X 10^ 




4 


0.30965978110985601234X 10^ 




5 


-0. 14256926805592749588 X 10^ 




5 


-0.49715203648917136888 xlO^ 




6 


0.16894039347875652311x10^ 




6 


0.63140021723358484451 X 10^ 




7 


-0.21057424914324066776x10^ 




7 


-0.90160716109902248718 x 10^ 




8 


0.2676311383571557235 x 10^ 




8 


0.96095356087432839558 X 10^ 




9 


-0.76173492536429366737X 10^ 




9 


-0.63690648044614848914 x 10^ 




10 


0.42006937485571 783327X 10^ 




10 


0.18830794813510467065x10^ 




11 


-0.12698208721358355433X 10^ 




11 


0.57163083968398816069x10^ 




12 


-0.53474107987358220271 x 10^ 




12 


0.19078924010583921336x10^ 




13 


-0.90835383726464555366X 10^ 




13 


-0. 1 1462564856718564988 x 10^ 




14 


-0.14158471202173284837X 10^ 




14 


-0.34295621170633088332x10^ 
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i 


J 




23 


15 


-0.83015840201854216866X 10^ 




16 


-0.11002212525564364623X 10'^ 




17 


-0.12227985275158903096X 10^ 




18 


-0.99733247054306510836X 10^ 




19 


-0.63431234698974736966 X 10^ 




20 


-0.26563714369463571347X 10^ 




21 


-0.64165338807344436646 X 10^ 


24 





0.60509199550873199769x10^ 




1 


-0.36606232348859748527X 10^ 




2 


0.12139375556095224965x10^ 




3 


-0.29330923818458009919X 10^ 




4 


0.55510496259245075635 x 10^ 




5 


-0.89523794399456470273 X 10^ 




6 


0.13690056986421268084X 10^ 




7 

8 


-0.14950913041347730905x10^ 
0.20276027262469305424x10^ 




9 


-0.24868958310934680048X 10^ 




10 


0.44658906685524416389 x 10^ 




11 


-0.50473239786133670714X 10^ 




12 


-0.20583143731838543317X 10^ 




13 


-0.31319042987920023734X 10^ 




14 


0.51903724071582528995 x 10^ 




15 


0.12406518345377091318x10^ 




16 


0.23558257104359627701 x lO'^ 




1 '~7 

17 


0.29049027840298512019x10 




18 


0.29832956761808455922 x 10^ 




19 


0.23110049545479607768 x 10^ 




20 


0.13877345160856433213x10^ 




21 


0.55381177327087428421 x lO'' 




22 


0.12541407130479812622x10^ 



4. Analysis of the series 

To obtain estimates of tlie critical properties, specially the critical line, from the 
supercritical series for the ultimate survival probability as given by the equation fl2T]) . 
we initially use d-log Pade approximants. These approximants are defined as ratios of 
two polynomials 
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Li 


iVi 




n 
U 


1 n 

iU 


1 n 
iU 


U.OU040 




1 1 
ii 


ii 


U.DUd4d 


n 1 


1 


1 






11 


11 


0.62266 


0.7 


10 


10 


0.85353 




11 


11 


0.84513 


0.8 


10 


10 


0.96256 




11 


11 


0.94246 



Table 2. Estimates for critical points obtained by d-log Pade approximants. Note 
that as the value oi D = D/{1 + D) grows the dispersion in the value estimates also 
increases. 

In our case the function /(A) represents the series for ^ In Poo (A). As /(A) is a function 
of one variable, we fix the value of D to calculate these approximants. For a fixed 
value of JD one pole of the approximant F will correspond to the critical point while 
the associated residue will be the critical exponent [3. We calculate approximants with 
L + M < 24, restricting our calculation to diagonal or close to diagonal approximants, 
which usually display a better convergence. Thus L = M + ^, where ^ = 0, ±1 and 
with D = D/{1 + D) ranging between and 0.8. Examples of estimates for the critical 
values of fi obtained from these approximants is given in Table [2] for different values of 
the diffusion. 

For each value of the diffusion rate, we calculate about eight approximants, 
obtaining the estimate of fic associated to diffusion as an arithmetic average of results 
furnished by these set of approximants and error bar associated to it corresponds to 
standard deviation of the estimates. From this we obtain the phase diagram shown 
in the figure [2l With the purpose of comparison with the results coming from the 
conservative simulations we use the variables a = fi/2 and -De// = C(D/{1 + aD). 
We also calculate the /3 exponent associated to the order parameter for differents values 
of D exhibited in the figure [3l Similar to what happens with the critical point a growing 
fluctuation for high diffusion rate values is visible. 

Turning back to the discussion about the critical line, we see that for higher values 
of the diffusion the dispersion increases, and estimates with larger error bars are found. 
Nevertheless, the exponent = 4 seems to describe well the calculated points of the 
critical line. However, this result would be different if we used approximants to series 
closer to the infinite diffusion rate limit. This error in the approximants for high values 
of the diffusion rate is attributed to the alternated sign of the series terms [15] . Another 
explanation would come from the fact that in the neighborhood of a multicritical point 
the reduction of a two-variable series to one variable leads to very poor estimates of 
the critical properties [12] close to a multicritical point. To overcome this problem we 
analyze the series using Partial Differential Approximants (PDA's) [T6|, that generalize 
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Figure 2. At left panel we have the phase diagram obtained by simulations (solid 
curve) and through Fades approximants for the supercritical series (circles). At right 
the log-log plot of the same quantity is plotted showing that value of the crossover 
exponent = 4 is a reasonable estimate for this data. 




Figure 3. Estimates of the exponent (3 for different diffusion rates. The dashed line 
indicates the DP value for this exponent. 



the d-log Pade approximants for a two-variable series. These approximants are defined 
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by the following equation 

Pl(x, y) = QMix, 1/)^^^ + ^n(x, 2/)^^^, (23) 

where P, Q, and R are polynomials in the variables x and y with the set of nonzero 
coefficients L, M, and N, respectively. The coefficients of the polynomials are obtained 
by substitution of the series expansion of the quantity which is going to be analyzed 

f{x,y)= ^ f{k,k')xY (24) 

k,k'=0 

into equation (l23l) and requiring the equality to hold for a set of indexes defined as 
K. This procedure leads to a system of linear equations for the coefficients of the 
polynomials, and since the coefficients fk,k' of the series are known for a finite set of 
indexes this places an upper limit to the number of coefficients in the polynomials. 
Since the number of equations has to match the number of unknown coefficients, the 
numbers of elements in each set must satisfy K = L-\-M-\-N — 1 (one coefficient is 
fixed arbitrarily). An additional issue, which is not present in the one- variable case, is 
the symmetry of the polynomials. Two frequently used options are the triangular and 
the rectangular arrays of coefficients. The choice of these symmetries may be related 
to the symmetry of the series itself [H]. In our case, the series presents a triangular 
symmetry when express in terms of the variables x = a and y = aD/ (1 + aD). 

It is possible to show [H] that we can determine the multicritical properties using 
the equation [23] and the hypothesis of that in the neighborhood of the multicritical 
point, the function / obeys the following scaling form 

f{x,y)^\Axrz(}^], (25) 



where 



and 



\Ax\ 



Ax = {x - Xc) - {y - yc)/e2, (26) 



= {y - Vc) - ei{x - Xc). (27) 

Here u is the critical exponent of the quantity described by / when Ay = 0, ei and 62 
are the scaling slopes [TB] and is the crossover exponent. 

On the other side, our calculation was successful when we use the method of the 
characteristics to integrate equation (!23l) . This is made by introducing a timelike variable 
r, so that a family of curves is obtained in the plane (x(r),?/(r)) (the characteristics). 
Such curves obey to the equations 

^ = QA4{x{r),y{T)), 
dr 

^ = RN{x{r),y{r)). (28) 

It is possible to show that integrating the equations (128!) from a specific point of the 
critical line, the resulting characteristic is equivalent to the the critical line. In figure 
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m we show a comparison between a characteristic curve and the simulational result 
[TT] . The number of elements in the sets of the calculated approximants was varying as 
follows: K = 55 - 190, M = = 20 - 53 and L = 15 - 54. 

In each of these curves, we calculate his inclination in the neighborhood of the 
multicritical point, determining a value for the exponent and the mean value for this 
exponent results as = 4.02±0.13, consistent with simulations in the particle conserving 
ensemble |llj . 




Figure 4. Comparison between a characteristic curve (solid line) and a numerical 
simulation result (circles). The coincidence is evident, including the region close to the 
infinite diffusion limit (inset). 

Using all the characteristic curves calculated we derive an 'average curve', 
calculating for each value of a on arithmetic average for -De//- This curve is shown 
in the figure [5] jointly with the result originating from the simulation and with the 
scaling form (1 — -De//) ~ (1 — a)^. In the same figure we see that the exponent = 4 is 
well fitted to the results of simulation and of the PDA's. This scaling form is based on 
the argument of the scaling function Z presented in the equation ( l25l) . where = 4 and 
Zq is a parameter properly chosen. We remark that this scaling form coincides with the 
characteristic curve and with the simulation even in the weak diffusion regime. This is 
somewhat surprising since its validity would be expected only in the neighborhood of 
the multicritical point. 

Unfortunately, even using the algorithm proposed by Styer p!8] we were not able 
to obtain precise estimates for the crossover exponent from the scaling form shown 
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Figure 5. Comparison between a cliaracteristic 'average curve' (solid line), the 
numerical simulation results (triangles) and the curve obtained from the scaling form 
(1 — Deff) ~ (1 — a)^(circles). At the right panel we see that = 4 seems to be a 
good estimative for the crossover exponent. 



in equation (125!) . However, integrating a set of approximants, we could determine the 
characteristic curves whose initial point is coincident with the critical point of the CP 
without diffusion of particles. These curves are estimates for the critical line of the CP 
model with diffusion going beyond the values achieved in [15] and [H] and corroborating 
the initial result of this last reference in that (p ^ 4. 
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5. Exact solution for finite systems 

An interesting approacli to study stocliastic systems is tlie exact solution of models 
with increasing numbers of sites L followed by extrapolations to the thermodynamic 
limit [5]. This is accomplished by considering operator S and its eigenvalues. For 
process governed by the master equation vanishing eigenvalues, /iq = 0, of the evolution 
operator correspond to stationary states of the system. For finite systems with absorbing 
states, only these states are stationary, and no active stationary state is found. To study 
the transition between an active and a stationary state, we may consider the behavior 
of the eigenvalue with the second smallest absolute value F = \ni\ of the operator. 
This eigenvalue is related to the quasi- stationary state [21] and will eventually become 
degenerate with /xq iii the thermodynamic limit, originating the phase transition. 

In a one-dimensional lattice with L sites, we construct a set of configurations that 
works as basis of the operator. Using periodic boundary conditions some of these 
configurations will be related for a symmetry Cl over cyclic rotations. These symmetries 
reduces the number of the independent vectors of the basis. 

On the other hand, a scaling behavior for F, with D held fixed, is given by the 
expression 



we may estimate the critical point ac{L) finding the intersection of the curves Yl and 
Yl+i [Sj. This procedure resembles the phenomenological renormalization group |22j . 
The sequence a^L) of estimates for a given value of D, with L ranging between 4 and 
14, was extrapolated to the thermodynamic limit, producing a critical line shown in the 
figure [6l Although in the higher diffusion region the coincidence between the two lines, 
shown in the figure, reduces there seems to be no doubt that = 4 is a good estimate 
for the crossover exponent. Actually, the estimative for this exponent using the exact 
diagonalization approach is = 3.87 ± 0.17. Probably for high diffusions is necessary 
to study larger lattices to determine more precise results for the critical point as well 
as its exponent. The exponent z is shown as function of the diffusion rate in the figure 
[71 For smaller diffusion rate, the exponent value coincide with that predicted by the 
DP universality class. However, for D > 0.2 occurs a growing of this estimative and a 
simple quadratic extrapolation gives z = 1.93, near to multicritical point, which is in 
reasonable agreement with mean- field value for this exponent, z = 2. 




(29) 



where A = |ac — a;|. Defining the quantity 




(30) 
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In(l-a) 



Figure 6. Comparison between a 'average' cliaraeteristic curve (solid line) and the 
exact diagonalization result (squares). At the right panel we see that </) = 4 is a good 
estimative for the crossover exponent. 

6. Conclusion 

We study the difussive contact process in one-dimension analyzing its critical properties. 
A special attention is devoted for the determination of the crossover between a classical 
mean-field behavior in the infinite diffusion limit and DP universality class presents when 
the diffusion rate is finite. To accomplish this task we use supercritical pertubative series 
and exact solutions for finite-size lattices. 

Calculating a supercritical series for the ultimate survival probability and analyzing 
it using PDA's we obtain estimates for the critical line of the CP model with diffusion. 
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Figure 7. Estimative for tlie z exponent calculated by some values of the diffusion 
rate. The dashed line indicates the DP value. 



The critical line was derived through the integration of the equation ( l23i) by the method 
of the characteristics. Direct results for the value of the crossover exponent using the 
scaling form [25] using Styer's algorithm [18] we were not be able to be obtain with 
an acceptable precision. However the method of the characteristics permitted us to 
calculate the critical line and, consequently, the value for the crossover exponent 0. Our 
result, = 4.02 ± 0.13 is in agreement with that derived in [IT] and explores a region 
of diffusion very close to the multicritical point. 

The technique of the two- variable supercritical series associated with PDA analysis 
was shown to be accurate enough to determine the critical properties in similar models 
|12] . Therefore, we believe that a natural extension for this work is analyze related 
models that apparently possess non-trivial multicritical points. This seems to be the 
cases of the pair-creation and triplet-creation models with diffusion, also studied in |llj . 
This research is already in course. 

Finally, the results provided by the exact diagonalization of the time evolution 
operator are according to simulational data as well as that ones obtained by the series 
analysis. Actually, in higher diffusion region, a greater fluctuation take place, suggesting 
that this system must be studied in larger sizes. Nevertheless, the critical curve indicates 
that ^ 4 is a good estimate for the crossover exponent. 
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